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ABSTRACT 


An algorithm, based on a series of analyses by Carlson, is presented for 
the symbolic integration of a class of algebraic functions. This class consists 
of functions made up of rational expressions of an integration variable x and 
square roots of polynomials, trigonometric and hyperbolic functions of x. The 
algorithm is constituted of four major components, viz., (i) reduction of input 
integrand to canonical form, ( ii) intermediate internal representation of 
integral, (iii) classification of output, and (iv) reduction and simplification 
of output to well-known functions. In the oral presentation, the algorithmic 
outline as well as some simple examples will be described. 
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INTRODUCTION 


Symbolic integration is concerned with the computer automation of 
analytic or closed form integration. Recently Moses (Ref. 6) presented an 
excellent survey of the last decade of efforts in this area, from which he 
discerned three main streams of interests, viz., from artificial intelligence, 
from algebraic manipulation and from mathematics. For our discussion here we 
shall characterize these interests oy two approaches, viz., the heuristic and 
the algorithmic. The former approach involves basically the application, by 
educated guess, of a class of methods each applicable to a specialized class 
of problems. The latter approach involves in the study of a general class 
of functions (e.g., rational functions) and attempts to determine if a given 
member of this class can be integrated in terms of a member of the same class 
or a well-defined extended class (e.g., rational functions extended oy 
logarithms). On the theoretical side, Risch's result (Ref, 9) stands out as 
commanding paramount importance, both in its generality and rigor. On the 
practical side, Moses' implementation stands out as the most successful to 
date. Therefore we shall describe our proposed approach in relation to the 
work of both Risch and Moses, emphasizing the parts which our approach will 
yield that are not contained in that of Risch and Moses, To this end we 
first highlight the contributions of these two investigators. 

Risch's algorithm is based on his reformulation of Liouville's theorem. 

It involves in first examining the hierarchical structure of the integrand, 
i.e., the finite number of elementary extensions that build it up from a 
ground field of rational functions. Once such hierarchy has been determined, 
the integration problem reduces to a systematic comparison of coefficients 
through differentiation. For our purpose here it is convenient to refer to 
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the integrand and the resultant integral as belonging to some input and 
output classes of functions. Using such terminology we should make clear 
that Risch's algorithm deals with only input and output that are elementary 
functions, that is, any of the following four clasttes, viz., (i) algebraic; 
(il) exponential, trigonometric and hyperbolic (as rational functions of 
exponentials); (iii) logarithmic; (iv) inverse trigonometric and inverse 
hyperbolic (as logarithms of algebraic functions). Moreover, the current 
state of Risch's algorithm handles only inputs that are rational expressions 
involving cases (li) and (iii). This point is made very clear in his most 
recent publication (Ref, 10). It is well known that many simple combina- 
tions of elementary functions do not possess elementary integrals, e.g., 

2 

X X 

e , e /x, etc. In such cases it is necessary to extend Risch's algorithm 
so as to enlarge either the input or the output class, or both. 

Another important point to note is that Risch's algorithm is designed 
mainly for the problem of indefinite integration. It is well known that 
many expressions are only integrable for a specialized interval, e.g., 
[0,2tt], ro,«], etc. Such definite integrations that possess no indefinite 
counterparts fall outside the realm of Risch's algorithm. Wang's (Ref. 11) 
investigations represent pioneering efforts in such definite symbolic 
integrations. 

Turning to the practical side of symbolic integration, Moses' program 
SIN is certainly the most general and sophisticated to date and is in our 
opinion the only viable one. SIN is, roughly speaking, constituted of 
three components: one made up of a collection of specialized heuristic 
methods, one consisting of the capability for rational function integration, 
and one implementing Risch's results. 
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With the foregoing preliminaries, we are ready to describe an algorithmic 

approach to handle a class of algebraic integrands. We shall label this 

approach Carlson's algorithm because it is based on a series of analyses 

rendered by Carlson and his associates in the last ten years (Refs. 2, 3> S 8> 

2 

and 12). The class of integrands is of the form r(x,y), where y is a 
polynomial in x, and r a rational function in x and y. This is the type of 
integrand that classically led to the study of elliptic integrals. At first 
glance this is a rather restricted class of algebraic functions. But in fact 
many trigonometric and hyperbolic integrands reduce to this form. The 
richness of this class of integrands is exemplified by a recently published 
handbook of 3000 integral formulas (Ref, 1), Our proposed approach will 
cover fifty to seventy percent of the items in the handbook. Furthermore 
the non-classical approach we shall describe holds great promise of 
developing to the case where definite integrals can be evaluated in terms of 
a host of other well-known functions (e.g. , Bessel and Legendre). It is 
important to make clear what we aspire to achieve that is not possible with 
Risch's algorithm. For this purpose let us divide our approach into four 
cases: (i) indefinite integration - elementary output; (ii) indefinite 

integration - non-elementary output; ( iii) definite integration - elementary 
output; ( iv) definite integration - non-elementary output, in all cases the 
input being algebraic functions. By definite integrations we mean those 
that possess no indefinite counterparts. It is clear from the previous 
discussions that at its current stage, Risch's algorithm cannot handle any 
of these cases, and even in its prospective completion it will only handle 
our case (i) here. On the other hand Risch's algorithm is far more general 
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in scope than our proposed one. Therefore the latter should be considered 
an adjunction or supplement, albeit an important one, to the general 
capability in symbolic integration. 

It is appropriate here to give some examples of the types of integrals 
that fall in the class under investigation: 

( 1 ) 

( 2 ) 

( 3 ) 

(*0 


Jt‘^(l+t 3 )'^dt 

Jt 2 (t 1< +at 2 +b)" 3/,2 dt 

o p _ l 

j ( A+B cos0 + C sinG + D cos 9 + E sin 9 + F sinP cos P)" ? dP 

o 

n (A sinh 9 ± B sinhP + C) do 

(E sinhG + F) 2 (G sinhG + H coshG)? 


It is also appropriate here to delimit the scope of our output. We are 
interested in adding to the output only three new functions to the elementary 
functions. These functions play the same role as the logarithm in the inte- 
gration of rational functions. These functions can be chosen as Legendre’s 
well known canonical integrals F, E, and rr or Carlson's standard functions 
(cf. Kef. 12). 

CARLSON'S ALGORITHM 

2 

The class of integrands, viz., r(x,y), where y = P n (x), has been the 
subject of intense study for almost two centuries. It is the simplest class 
beyond rational functions, but even such a class frequently leads to non- 
elementary integrals. For n = 0, 1, or 2, the integral is elementary; for 
n = 3 or U, the integral is elementary or elliptic in which case up to 
three new functions need be introduced; and for n ? 5, the integral is called 
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hyperelliptic, in which case more than three new functions are required. Two 
main classical approaches to handle this class were due to Legendre and 
JacoDi, respectively (Refs. 1, F), and a more recent one due to Carlson and 
his associates. We believe that Carlson's approach is superior to the other 
two for the purpose of symbolic integration. Such belief is based on the 
observation that Carlson's representation possesses certain symmetry that 
offers significant simplification in formal manipulation and economy in 
storage of formulae. Here we give a few examples to illustrate this remark, 
reserving for a future report a detailed comparison of the three approaches. 
First, many integrals involving different Jacobian elliptic functions, say, 
cn u, nc u, dn u, nd u, cd u, dc u appear in Carlson's representation as 
various permutations of arguments of the same integral. Thus in such a 
situation one reduction formula is required for Carlson's case as compared 
to six in Jacobi's case, (see remark and example in Ref. 8 , p. 327). Second, 
in the case of reduction of Legendre’s integrals, different substitutions 
are needed according to various types of roots of P n (x). Here one either 
heuristically chooses the substitution (which is not always successful, see 
Ref. 7, Pp. 600-602), or stores rather substantial tables of transformation 
formulae (see Ref, 13, Pp. 78-79, 8 ^- 91 ), in addition to the regular 
reduction formulae. With Carlson's approach this procedure is not necessary 
because of the possibility of rendering the integrand into a canonical form. 
Third, a number of linear and quadratic transformation formulae are needed 
in Legendre's and Jacobi's approaches (see Ref. 1, Pp. 12-lb, 38 -Ul). In 
the present context these formulae would be required for algebraic 
simplification, which is crucial to the success of any program in symbolic 
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computation. With Carlson's approach the linear transformations are rendered 
trivial by permutation symmetry while the quadratic transformations take on 
a canonical form (cf. Ref. 12, Pp. 203~208) . 

In this section we shall translate Carlson's results into an algorithmic 
scheme. This scheme can be summarized by the following steps: (l) Reduce the 
input function to a canonical form; ( 2 ) represent the output integral formally 
as a certain Carlson’s R-f unction (which is just an intermediate representa- 
tion); ( 3 ) examine the parameters of the R-f unction to characterize the out- 
put, for example, whether it is elementary or elliptic; (U) reduce the 
output to a combination of elementary functions and well-known elliptic 
integrals. 

Here we are mainly interested in the "algorithmic outline" and there- 
fore to avoid substantial digressions into mathematical details, we reserve 
the description of Carlson's R-function to Appendix I. It suffices here to 
specify the notation R(a:b 1 , . . . .b^jz^, . . . ,z^) to indicate that in general 
this function is characterized by k variables and k+1 parameters. For the 
major bulk of our work here k is between 2 and t. For brevity we shall 
refer to this function as R and indicate whatever change in parameter or 
variable oy specifying that parameter only, e.g., R(b>l) s R(a; b^, . , . , b^+1, 


Procedure I: Canonical Reduction of Input 

1,1 Given input expression r(x,y); if it is already in the form of 

linear factors r(x,y) = fj (g^x+w^), go directly to Procedure II, 


i=l 


else proceed. 
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I.c' Reduce to linear fraction r(x,y) -♦ [P 1 (x)+yP ? (x) V[P^(x)+yP l4 (x) ], 
where the P^(x) are polynomials in xj this reduction is possible 
because in any polynomial in x and y, the even powers of y can be 
expressed as straight polynomials of x and the odd powers as 
polynomials of x multiplied by y. 

1.3 Multiply both the numerator and denominator by T P^(x)-yP lj {x) ] and 

then by y to transform it to the form r 1 (x)-t-r ? (x)/y , where r^ and 
r ^ are i Honal function? of x. 

l.U Integrate r^(x) by standard rational function integrator, 

1.5 Apply square-free and partial fraction decomposition to r„(x) , 

A A 

resulting in r^x) =^-4-...+ -* (see, e,g., Ref. 6, p. 551)* where 

1 s t 

the S.'s have only simple roots and are relatively prime to each 

V 

other. 

1.6 Factori-e each S into linear form 77 (x+z ), wtere z. may be real 

J= 1 J <J 

or complex. Now the integration of r 2 (x)/y reduces to that of 

Aj/(S^y) which in turn reduces to that of x ra y 1 77 (x+z^) f . There 

J —1 

is of course the practical problem of factorization. But the 
problem here is no worse than the case for rational function 
integration. Furthermore, for most cases of practical interest 
the S,'s are polynomials of very low degrees, as evidenced by the 
several thousand formulae in Ref. 1. 

Procedure II: Intermediate Representation of Integral 

II. 1 lor the integration of algebraic functions, extreme care must be 
exerted to avoid the indeterminacy of integration across a branch 
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point, tor example, in the Legendre formulation, most of tne above- 
mentioned substitution* are required to avoid such orabigui v , In 
Carlson's formulation the situation is simpler. Here we need only 
consider the indefinite integration variable as either always to the 
left of the smallest root of y on the real axis or to the right of 
the largest real root, Bor definite integration we need to allow 
the degenerate cases of one or coth integration limits bein'?, branch 

points, for the sake of such distinction we need first to factorize 
k 

y -* n (x+zj. 

i=l 

11 . 2 Bor indefinite integration, a user of this algorithm is asked to 
indicate the appropriate integration range, i.e., whether it is to 
the left or the right. We then proceed to give either one of the 
following representations: 

k -b. ^ ^ 

ef“ /7(t+z.) x dt = R(a;£;x+z), (s) 

x i-1 1 

k -h. k -b, 

J* (7 (t+z i ) <it = x /7 (x+z A ) i p(l;a,^;l,lf |) , (6) 

0 i=l i*l 

where a ^ b l' * * in a = b l + ’ * * +b k~ 2 ’ in (6); t? = (b 1> ...,b k ), 

- (x+Zi.x+z^^.^x+z^, etc. 

11. 3 For definite integration, check for degeneracy; if none apply 2 
above, else use the following: 
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( 7 ) 


r x l. .a- 1 , xfl'-l k , “ b i 
c (t-x r (x.-t) fj (t+z.) dt 

X n i=2 


- (x -x r(a)r(« t ) 

' 1 0 ; r(a+a ' ) 


It -b , / ^ 

/7(x+zJ R(a;b;l, 
i»2 0 1 \ 



where a+a ' = b.,+***+b, . 

1 k 

II. k The R-function so obtained has a very special form, viz., the b^'s 
(and therefore a) are integers or half- integers. Since y is the 
square root of an n tl1 degree polynomial, there must be n parameters 
among a and b^'s that are ha If- integers, where h « n or n+1 for n 
even or odd, respectively. This is the so-called R-function of index 
h, which can always be expressed in terns of R-functiona of h 
variables in which a, b^, b 2*’*' ,i:i n-l are half-integers and b^ is 
an integer. The next step in the algorithm is to actually achieve 
■•■■hat via recurrences. First rearrange the R-function to the form 
such that are integers in increasing order and L = -m; 

this can always be done due to permutation symmetry; in most 
cases of practical interest, b^...b k should be small integers, say 
+1 and +2, 

II. 5 Recall from Step 1.6 that we now have integrals of a polynomial 
divided by S^Jy, i.e., <y m jx ra /(sjy)dx + w m . 1 Jx n ’" 1 /( S^y)dx + ..., 
where each m, m-1,... appear as the index b^ in the corresponding 
R-function. Recursive operation n-times with the identity 
cR = a'R(b^+l) + a z^ R(a+l,b k +l) (9) 

will combine these m terms and reduce each b to zero. Here z. 

K X 
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II. 6 


II. 7 


is the k tn variable jf the R-function znd the identity is readily 

derived from the general recurrence relation, eq, A(ll), 

Now b. ,..,b, , ari all exactly t , which makes it efficient for the 
h+l k-l 

following recursive operations 

(o-l)(z K _ 1 -z k _ 2 )R - (c-DrRCa-l^^-D-RCa-l.b^g-Dl (9) 

Here z. ^ z 0 because S has only simple roots. 

K*1 K"c jt 

The above two steps reduce the integral CA^/(S^y)dx to an R-function 
of index h, i.e. , 

R(a;b lt . . . , b^ ,0 , , * * * 0 J » * * . . z^) R{ a ^ , * * * * ^h* ^1 * * * 9 * ^h^ * ( 1^) 


Procedure III: Classification of Out-nt 

The classification depends on the parameter h -2, h f 6... 

111.1 If h=0, that means the input is a rational function of x only and we 
have only come this far by mistake. Here an error return with such 
a message would be appropriate, 

111. 2 If h=2, the integral is elementary and each R-function is expressible 
in terras of two "basis" functions, viz., 

and (ii) 

_ i . . - i 

( 12 ) 


R( 5 ; 2 , 1 ;S 1 ,z 2 ) * (z 2 -z 1 )"^cos*‘ 1 (z 1 /z 2 )^. 


This case would have been handled more easily by rational substi- 
tutions, but here we also want to "catch" integrals that have a 
misleading non-elementary appearance, and so this capability is 
retained here. Now set -flag P = true and apply Procedure IV. 2. 
(Note: The flag is to distinguish between the pseudoelliptic and 

elliptic cases.) 
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III. 3 The cose h=h constitutes the major bulk of our effort here. For example 
at least two thousand formulae in Ref, 1 belong in this category. 

(i) In general a function of this class is of the form Rfaj^.b^, 

(13) 

( 1*0 

(15) 

(16) 

symmetric in the variables z^, z^, and Zy 

(ii) There are two important special cases caused by one or two of 

the following conditions: (a) any z^ = 0, (b) b^ = 0, and (c) any z^ 

= z., Here a flag k * U should be set and for any occurrence of the 
J 

above conditions, decrease k by 1 and reduce the R-function by one 
of equations (2b) - (27). Now go to the three cases according to k = 


b and is reducible to h basis functions, viz,, 

(*1*2*3)^ S 

VVVV s VV’ 

H G^l’W S r ("z5^i^2*z 1 »\,»^3)* 

VV*3»Zl,)' 

The last three functions are analogous to Legendre's canonical 
integrals F, E and tr (Ref. 1 , p. 8), All these functions are 


A A 


2, 3, *». 

(iii) Special case k=2: R has been reduced to R(a ;b^,bg;z^,z^) . 


(a) If b^ or b^ is integer, we have pseudoelliptic case. Now 
set P = true and apply Procedure IV. 2. 

(b) Else the answer is expressible in the form CR = Cj^ + 
where the C's are polynomials in z^ and z^, and 

k k * r e 5 (17) 
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Ill.b 


They are analogous to Legendre's complete integrals K and K. 

Now set P 3 false and apply Procedure IV, t? to obtain this answer, 
(iv) Special case k— 3 i R has been reduced to R(a;b^,b^,b^;£^ t £ l1 ^) . 
The answer is expressible in the form 

( 18 ) 


cs . Cft. > c 0 B a . C A (£,S^)^. 


where the C's are polynomials in z, , 'z ^ and z . Apply Procedure 

J. C j 

IV. 3 to obtain this answer. 

(v) General case K=A, The answer is expressible in the form 

CR * W * Vo * Vh + C *^* J 3» 4 - < J 9) 

A A A 

where the C*s are polynomials in z . , and z , Apply Procedure 

■L c j 

IV. to obtain this answer. 

The cases h = o, 8, ... in general deal with hyperelliptic integrals. 
In practice there are many such integrals that hove hyperelliptic 
appearance but are actually reducible to the elliptic case (e.g. , Kef. 
1, Pp. ^Sl-c’71) . We are only interested in such degenerate cases. As 
before first set a flag k=h. 

(i) First check lor obvious degeneracies (a) = 0, (b) - 0, or 

(c) z^ = Zy Each occurrence of any of these conditions reduces 
k by 1. 

(ii) For each equal pair (z.,z.) compute the new parameter b. . = t>, 

i J ij 

+■ bj. Collect all integer b^’s ar| d now apply the some ooere- 
tions as in II. k - II. 7 to further reduce the R- function. 

The reduction of each parameter is indexed by a corresponding 
decrease of k by 1. 

(iii) Now if k < S, go to step m. 3» else return the hyperelliptic 11- 
function with no further reduction. A typical user may still 
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find this R-function to be useful because many properties are known 
about it (Ref. ?> 


Procedure IV. Reduction of R-function to tasis Functions 

This procedure consists of three sub-procedures (k<?, 3» M for the reduction 

of R(a;b ,b , , . . ,b ;z , . . . ,z ) to a set of basis functions. The basic technique 
1 <. K 1 K 

in this Procedure involves the application of the recurrence formulae (;_•’»)- (73 ) . 
tor this preliminary report we shall only present the case k=P, leaving the 
more tedious ones for the final detailed report. 


IV . 2 Reduction of 7-parameter R-functions: 

If the flag P - true, the two basis functions are elementary (ll1-(lP}, 
else they are R^. and R^ (T7V*(lB} . For this Procedure we need the 
following four formulas which may be derived from (?Q)-(3t )• 
(c-l)R(b.-l) - (a'-l)R * a^Rta+l); 


(bO) 

(Tl) 

(-) 


^(c-ll^bj-l) - z j (c-l)R(b i -l) - (b J -z i )(o , -l)R; 

az.z.R(a+l) = a'R(a-l) - (b.-a)z.R - (b.-a)z.R 
i j i i J J 

a ' b . (z -z . )R(b , +l) = c {(b -a) ( z ,-z )R+z (c-l)f R-R(b . -1) 1; 

Xljl 1 J ^ J 1 

Now the algorithmic steps are as follows: 

(i) Given an R-function R(tj;b^;b and the busis functions 
R^ and R^; 

(ii) apply eq. (p?) repeatedly to R^ and R„, to obtain 

R(a;b i ,b 2 ;zi > z 2 ) and R^+ijb^b^z.^) ; (bM 

( iii) apply eq. (pf>) to compute R(a -.b^-l.b.^z^z^) ; 

(iv) recur on b^, with eq. (?3) to compute R(a ; b 1 ~l ;b •,z 1 ,z. 1 ) ; 

(v) apply eq. (?'l) to compute R(a » b 2 " 1;z i> z ^ ; 
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(vi) recur on b_, with eq, (?3) to obtain R{a;b. ,b„;z. ,z,,) ; 

c X c X 4- 

(vii) return result in the form C^R^ + C 2 R 2 * 

Procedure V, Transformation of Output to Legendre's Notation 
This Procedure is used only optionally for those who prefer Legendre's 
notations. The implementation of such transformation is basically clerical 
in nature, with the appropriate formula given by Carlson (Ref, 3 ), There- 
fore no algorithmic flow is indicated here. 
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APPENDIX I Properties of R-functions 

The R-function is essentially a generalized nypcrgeoraetric function of n- 
variaoles. It is symmetric with respect to simultaneous permutation of the 
u-parameters and z-variables. For detail properties see Carlson (Ref. ? ). 
Here we summarize those properties useful in the text. 

1 . Special values of parameters and variables; 


R( a ; » , , , » i » . . . * z^) = R(® ?b^, , , , ,b^_^}z^, , . . ,z _-) 


B(a,a’-b ) 


R(a;b 1 ,...,b tj ;z 1 ,...,z n _ r O) - ^ a ; a i 


R(a;b 1# . . . ,b ; 


z -\ * • * * > z „-l ) t 


R(a ;b^, . . . ,b^ ; z^, , . , , ^n-l^ ~ R(o ib^, . . . t b n _^+b^ ; z^, , . . , 2 n _j ) > 


C = b^+'-'+b^, a* = C-a, and B(a,a') = r(a)T(a ' )/r(a+a ' ) , 


(Zb) 

(20) 

( 27 ) 

(26) 


Owing to permutation symmetry, the zero or equal pair of variables in the 
last three formulas is not confined to the last position. Moreover, 
these identities can obviously be generalized to multiple zeros or 
multiple pairs, or any P-tuplets of equal variables. 

2 . Recurrence relations: 


(c-l)R(b i ~l) = (a'-l)R + az i R(a+l), i=l,...,n (29) 

n 

cR « ][>.R(b + 1 ), (30) 

i=l 

cR(a-l) = z (31) 

i=l 

n 

acR(a+l) = C cR - a’R(b.+l)l. (32) 

i=l 111 
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